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A strong  need  exists  in  active  search  and  6urveiliar*ce  systems  for 
simultaneous  multichannel  signal  processing  capability  The  data  acquisition 
rate  and  its  processing  in  real  time  have  already  strained  the  capacity  of 
conventional  search  systems,  and  in  some  cases  an  answer  if  being  sought 
in  the  duplication  and  parallel  operation  of  complex  equipments  This  trend, 
although  dictated  under  certain  circumstances,  is  undesirable  from  the 
standpoint  of  economic  and  reliability  factors  obviously,  but  also  it  suffers 
from  limited  capability  to  evolve  into  more  flexible  concepts  in  the  near 
future  Fortunately  a number  of  signal  processing  techniques  have  been 
proposed  in  recent  years,  and  some  of  them  have  been  implemented  with 


reasonable  success 


' J The  panchromatic  optical  correlator  is  one  attempt  at  satisfying  a 
particular  requirement  of  high  data  rate  receivers  matched  to  the  transmitted 
waveform  The  technique  is  presently  etill  in  x breadboard  stage  and  its 
applicability  is  not  fully  proven  as  yet,  but  it  offers  a reasonably  promising 
expectation  for  the  simultaneous  processing  of  variously  Doppler  distorted 
signals  This  simultaneous  approach  would  obviate  the  need  for  a sequential 
scan  of  all  the  Doppler  possibilities,  which  customarily  is  mulMplexed  with  a 
real-time  correlation  process  The  technique  is  particularly  applicable  in 
^'low  Q ^ wide  bandwidth  systems  in  which  Doppler  distortions  cannot  be 

approximated  by  simple  frequency  shifts  of  the  signal  components  Instead 



the  true  scale  factor  relationship  must  be  maintained  between  frequency 
components  of  the  transmitted  and  received  signals 


where 


v 


target  radial  spued  ti«um«d 
constant  for  signal  duration  T , 


c “ signal  propagation  speed. 

a r Doppler  distortion  factor, 

w - any  particular  radian  frequency 

of  the  transmitted  signal,  and 

w corresponding  radian  frequency 

of  the  received  signal 

An  ideal  correlation  receiver  would  cross  correlate  the  received 
signal  not  only  against  a replica  of  the  transmitted  waveform  but  against 
the  entire  range  of  all  possible  Doppler  distorted  versions  of  this  replica 
(as  dictated  by  a likely  range  of  v values  in  Eq  (.1) ) Such  a procedure 
becomes  prohibitively  expensive  when  waveforms  of  long  duration  T are 
used  to  obtain  increased  processing  gain  A safe  criterion  for  avoiding 
significant  degradation  in  the  correlation  process  due  to  timebase  contractions 
or  expansions  of  the  type  represented  by  Eq  (1)  is  the  one  which  states  that 
timebase  errors  between  the  received  signal  waveform  and  the  reference 
function  with  which  it  is  correlated  shall  not  accumulate  over  the  signal 


duration  by  more  than  one  quarter  period  of  the  highest  frequency  contributing 
significant  energy  to  the  power  spectrum  of  the  signal  For  large  T , the 
result  is  a fantastically  large  requirement  on  the  number  of  distinguishable 
Doppler  channels  that  must  be  processed  simultaneously  to  cover  all  possible 
target  motions  A coarser  spacing  of  Doppler  channels  leads  not  so  much  to 
a lower  resolution  on  target  velocity  determination  as  to  an  outright  failure 
to  detect  the  target  Therefore  the  entire  range  of  Doppler  distorted  signals 
must  be  scanned  with  a resolution  which  is  proportional  t*>  the  signal  duration 
T Electronic  techniques  can  be  used  to  time  compress  the  original  received 
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signals  and  scan  the  entire  Doppler  range  while  processing  the  received 
signals  in  real  time.  However,  such  techniques  are  not  generally  well 
adapted  to  multichannel  inputs,  as  from  many  receptors.  If  time -compression 
is  used  to  permit  a scan  of  Doppler  distortions,  it  is  not  then  available  to 
allow  a multiplexing  of  various  received  signals  onto  a single  timebase  so 
that  a single  electronic  correlator  can  perform  the  receiver -multiplexed 
correlation  Moreover,  the  number  of  received  signals  to  be  cross  correlated 
frequently  is  great  enough  to  require  multiple  correlation  channels  if  real- 
time operation  is  necessary  In  such  cases,  optical  correlators  have  an 
obvious  appeal  in  their  ability  to  correlate  a large  number  of  channels  simul- 
taneously without  duplication  of  equipment  The  panchromatic  optical  cor  - 
relator  takes  the  further  step  of  simultaneous  processing  all  Doppler  dis 
tortions  while  maintaining  the  multiple  channel  capability  at  the  input 

The  principles  of  coherent  optical  correlators  are  well  described  in 
the  literature  * Basically  an  optical  correlator  consists  of  one  of  several 
configurations  of  lenses  and  apertures  which  perform  the  following  generalised 
two-dimensional  averaging  process  (with  the  exception  of  the  Optical  Matched 
Filter,  shown  in  Fig.  2).- 

j(u  x + u>  y) 

H (w  . w ) «=■  JJ  f(x„y)  g(x.y)  e X y dxdy  (21 

y A 
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L J Cutrcna  E N Leith  C J Palermo  and  L J P^rcelloj 
"Optical  data  processing  and  filtering  systems,"  XRE,  Vol  IT  6 No  3 
(June  I960) 
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The  relation  among  coordinate*  of  Eq  (2)  ia  ihown  in  Fig  1.  and  the 

subscript  A indicates  integration  over  the  optical  aperture  The  weighted 

j(u>  x + w y) 

averaging,  in  which  the  kernel  e x ^ appears,  is  a property  of 
optical  system*  that  has  been  known  for  some  time  ^ Physically  the  kernel 
represents  the  phase  of  the  light  amplitude  vector  at  the  point  u>  ) 

determined  by  Huygen'a  principle  from  ths  light  amplitude  at  point  |x,y) 
and  the  integration  represents  a coherent  summing  of  light  amplitudes  over 
the  extent  of  one  of  the  aperture  planes  perpendicular  to  the  *»ptic  axis  in 
particular  various  restrictions  and  special  conditions  upon  Eq  (2)  uhich 
can  be  simulated  by  lenses  and  apertures  can  reduce  the  general  expression 
of  Eq  (2)  to  many  of  the  commonly  useful  transforms,  including  Fourier  and 
Laplace  transforms,  correlations,  convolutions,  power  spectra  ambiguity 
functions,  etc  In  addition  with  the  use  of  cylindrical  lenses  it  '#  possible 
to  reduce  the  two  dimensional  integration  of  Eq  (2)  to  s one  dimensional 
one,  in  which  case  the  unused  dimension  contains  a large  number  of  independent 
information  channels  which  can  be  processed  simultaneously 

A particularly  useful  configuration  for  performing  optical  correlations 
has  been  described  previously  and  is  shown  schematically  in  Fig  2 It 
uses  a coherent  monochromatic  light  source  to  project  the  Fourier  transform 
S(wx)  of  the  moving  signal  function  a(x]  x^)  upon  the  reference  function 
Fourier  transform,  R*(a>x)  , so  that  the  light  amplitude  distribution  trans- 

mitted through  the  frequency  plane  in  which  R*(*.  ^1  is  inserted  is  proportions! 


Max  Born.  Optik  (Julius  Springer.  Berlin.  19311,  Sec  46  4? 


to  S(u»x)  R*(^x)  An  inverse  Fourier  transform  performed  by  the  lens 
to  the  right  in  Fig  Z produces  the  desired  correlation  function  4><x  ) in 
the  output  plane  This  optical  arrangement  is  called  a matched  filter 
because  it  uses  in  analogy  to  electronic  matched  filters  a spatial  frequency 
filter  R*(u»x)  matched  to  the  signal  «(Xj)  In  this,  or  other  configu 
rations  » different  reference  function  must  be  inserted  for  each  different 
Doppler  distortion  of  the  received  signal  Therefore  to  search  the  entire 
range  of  possible  target  motions  the  reference  function  must  be  changed 

I continuously  and  the  various  correlation  functions  for  the  different  Dopplers 

performed  sequentially 

The  panchromatic  optical  correlator  depends  for  its  operation  upon 
the  fact  that  a Fourier  tranaform  relationship  exists  between  the  light  amplitude 
distributions  in  planes  located  with  respect  to  a lens  ar  shown  in  Fig  Z 
This  relationship  has  been  known  for  some  time,  but  its  practical  application 
to  optical  correlators  first  found  expression  In  the  conventional  and  matched 
filter  configurations  It  is  further  exploited  in  the  panchromatic  correlator 
where  advantage  is  taken  of  the  fact  that  the  scale  factor  connecting  the 
coordinate  of  a frequency  plane  and  the  corresponding  spatial  frequency 
depends  on  the  light  wavelength  ' 


where 


Zw 

“ • TT  x2 


(3) 


w = radian  spatial  frequency  associated  with  x, 
X - light  wavelength, 
f lens  focal  length  and 

x^  coordinate  defined  on  j.’tg  Z, 


It  is  obvious  from  Eq  (3V  by  its  dependence  upon  X that  s con- 
tinuous range  of  scale  factors  connecting  u>x  and  x^  can  be  obtained  by 
varying  the  color  of  the  light  source  In  particular  if  we  define  U>  » U>  ' viw 
and  relate  x ^ to  u>  Eq  (3)  can  be  made  to  simulate  Eq  (1)  wherein  the 
reflected  signal  frequency  components  appear  as  scale  factor  distortions  of 
the  original  signal  components  Thus  a slight  broadening  o f the  spectral 
width  of  the  light  source  color  filter  uhich  is  made  as  monochromatic  as 
possible  for  usual  correlator  configurations  can  simulate  in  the  panchromatic 
correlator  the  Doppler  distortions  imposed  upon  the  signal  bv  reflection  from 
a moving  target 

The  complete  configuration  of  a panchromatic  correlator  is  shown  in 
Fig  3 It  differs  from  earlier  correlator  types  by  including  a broadened 
spectral  width  In  the  light  source  and  filter,  a diffraction  grating  to  allow  the 
uae  of  a heterodyned  version  of  the  received  signal,  a single  fixed  reference 
function  inserted  in  a frequency  plane,  and  a color  separation  prism  to  resolve 
the  various  Doppler  channels  in  the  output  plane  The  reference  function,  or 
replica  of  the  transmitted  signal  spectrum,  is  introduced  as  a matched  filter, 
i e , the  complex  conjugate  of  the  signal  Fourier  transform  i*  inserted  in  a 
frequency  plane  as  a filter  which  matches,  for  a given  center  wavelength  \° 
of  the  light  source,  the  undistorted  transmitted  signal  Neglecting  the  ef  ect 
of  the  heterodyning  operation  performed  on  the  received  signal  in  order  to 
conserve  bandwidth  prior  to  insertion  into  the  correlator  (to  be  described  in 
the  next  paragraph),  a continuous  range  of  Doppler  distortions  can  be  accom 
modated  by  the  continuous  spectrum  of  light  source  frequencies  without 
changing  the  physical  format  of  the  reference  function  This  format  is  fixed 
in  terms  of  the  x£  coordinate  of  Eq  (3),  but  Eq  (3)  shows  that  a range  of 
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corresponding  ^ , and  therefore  a . values  is  available  because  of 

the  range  of  X values  present  In  this  context,  the  panchromatic  correlator 
can  be  considered  to  be  a superposition  of  many  coherent,  monochromatic 
correlators  with  each  light  source  frequency  corresponding  to  a particular 
Doppler  distortion  of  the  transmitted  signal  The  various  light  frequencies 
are  separated  to  resolve  Doppler  by  the  prism  just  to  the  left  of  the  output 
plane  in  Fig  3 

A panchromatic  correlator  could  be  built  according  to  the  foregoing 
description,  and  presumably  it  would  properly  correlate  unheterodyned 
Doppler  distorted  signals  as  claimed  However,  for  any  signal  function 
whose  spectrum  does  not  extend  down  to  *ero  cycles  it  would  make  inefficient 
use  of  the  number  of  resolvable  elements,  or  time  bandwidth  product,  in  the 
optical  aperture  Any  optical  system  with  apertures  on  the  optic  axis  is 
necessarily  a low  pass  filter,  Hence,  better  use  of  the  available  bandwidth 
can  be  made  by  heterodyning  the  signal  spectrum  down  toward  tero  cycles 
The  heterodyning  operation  should  not  transpose  the  positive  frequency 
spectrum  into  the  negative  frequency  region,  or  vice  versa,  however,  for  this 
would  lead  to  overlays  of  positive  and  negative  frequency  components  and 
distortion  of  the  signal.  An  expression  analogous  to  Eq*  $1  and  i ? but  also 
containing  the  heterodyning  operation  would  be  written 

wx  3 aw  ‘ “h  = ' If  X2  ' -h  ’ <4) 

where  ox  ~ radian  frequency  shift  ejected  by  heterodyning  it  is  apparem 
from  Eq  (4?  that  the  simple  scale  factor  relating  snd  x^, 

in  Eq  (3i  no  longer  holds  in  Eq  (4i  Thus  the  heterodynint  operation 
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appears  to  eliminate  the  possibility  of  simulating  Doppler  distortions  by  a 
spread  in  the  light  spectrum  However,  by  inserting  a diffraction  grating 
immediately  after  the  first  objective  lens  of  Fig  3,  it  is  possible  to  aim 
ulate  Eq  (4)  so  that  a spread  in  light  frequencies  again  simulates  Doppler 
distortions  when  illuminating  a heterodyned  received  signal  function,  while 
the  bending  ol  the  optic  axis  into  the  first  order  image  of  the  diffraction 
grating  produces  the  constant  frequency  shift  Bor  signals  hetero- 

dyned down  near,  but  not  through,  tero  cycles,  the  scale  factor  relation 
between  and  x^  which  depends  upon  X , and  the  bending  of  the 

optic  axis,  also  dependent  upon  X , just  compensate  to  give  a resultant 
which  is  independent  of  X This  is  easily  shown  as  follows 

Figure  4 shows  the  bending  of  the  optic  axis  into  the  diffraction  orders 

3 

n = 0,  ± 1 The  usual  diffraction  equation  ' gives  the  locations  of  principal 
maxima  as 


where 


nX  = d sin  9 , t5) 

n = the  diffraction  order,  and 
d - slit  separation  distance  on  grating 


Also,  from  Fig  4, 

= 8 sin  e\ = 


(6) 


F A Jenkins  and  H E White,  Fundamentals  of  Physical  Optics 
(McGraw  Hill  Book  Co  , Inc  , New  Vork,  14.3*7)  1st  ed  , C"hap 
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using  Eq  (5)  Since  the  coordinates  of  Fig  3 are  measured  with  respect 
to  the  optic  axis  before  bending,  it  is  necessary  to  rewrite  Eq  (3)  in  terms 
of  distances  measured  from  the  optic  axi6  after  bending  It  is  shown  in 
Appendix  I that  the  relation  of  Eq  (3)  still  holds  in  the  first  order  diffracted 
image  provided  coordinates  are  measured  from  the  diffracted  optic  axis 
Therefore,  in  terms  of  the  coordinates  on  Fig  4, 

“*  * °“  - "h  ’ • If  [x2  - Vm]  • 0) 

where  x^(X)  i»  defined  on  Fig  4 
Substituting  from  Eq  (6;, 


Qu)  - IX).  - 


Zv 

~kT  X2 


Zvna 

fd 


(8) 


For  n = - 1 in  Eq.  (8),  the  term  represents  the  constant  u)^  , 

and  therefore  the  frequency  shift  effected  by  heterodyning  Then  th*»  ’•rmainder 
of  Eq.  (8)  becomes 


aw 


2n 

~e~  x- 
Kt  t 


or 


aXu  - 


f 


,o 

K It) 


(9) 


where  X°  is  that  wavelength  corresponding  to  a = 1 (no  Doppler  distortion) 
It  is  now  obvious  from  Eq  (9)  that  the  reference  function  can  be  constructed 
as  a function  of  x^  , which  is  related  to  the  transmit. ed  (undistorted)  signal 
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frequencies  a by  the  constant  - - x°  However,  the  reference 

function  will  also  represent  any  other  value  of  c»  for  which 


a X X 


The  ability  or  the  panchromatic  correlator  simultaneously  to  represent 
all  Doppler  distortions  of  the  transmittei  signal  (all  a ie  illustrated  graphically  in 
Fig  5 Here  the  Doppler  distorted  received  signal  au>  is  plotted  horizontally, 

X and  the  Doppler  distortion  factor  a are  plotted  vertically  The  frequency 
axis  for  the  transmitted  spectrum  <*>  is  plotted  horizontally  at  \ X°  and 
covers  the  bandwidth  from  A to  B For  a single  transmit  ed  component  w 
Eq  (10)  gives  a hyperbolic  relation  between  a and  K , which  xs  shown  in 
Fig  5 by  the  hyperbolic  curve  connecting  A , A and  A , all  corresponding 
to  minimum  (or  B’  , B and  B for  maximum  u>  Thus,  the  horizontal 
segments  A B A B , and  A B represent  the  received  signal  spectrum 
width  after  negative  zero  and  positive  Doppler  distortions  respectively  The 
reference  function  recorded  as  a function  of  for  the  transmitted  signal  of 

spectrum  w *th  A B accomm  dates  all  distortions  of  this  signal  represented 
by  other  horizontal  segments  connecting  the  two  hyperbolic  curies  of  Fig  ^ 
so  long  as  Eq  (10)  holds  true  The  light  source  spec'r  stn  w idth  is  cho«en  to 
insure  the  latter  The  constant  heterodyning  shift  -£lL8  appears  on 

Fig  5 as  a straigh*  vertical  line,  independent  of  X 

It  is  apparent  from  Fig  4 and  Eq  (8)  that  negative  x^,  and  n values 
simulate  the  Doppler  distortion  and  heterodyning  operations  that  have  occurred 
on  the  positive  u)  part  of  the  signal  spectrum,  while  positive  x and  n 
values  do  the  same  lor  the  negative  spectrum  Now  it  is  obvious  that  a mitchtd 
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filter  operation,  of  the  type  shown  in  Fig  2,  can  be  carried  out  by  inserting 
the  complex  conjugate  of  the  transmitted  signal  positive  frequency  spectrum 
in  the  ux  plane,  centered  a distance  fr&rn  the  original  axis  shown 

in  Fig  4,  and  the  complex  conjugate  of  the  negative  spectrum  symmetrically 
displaced  a distance  ) from  the  original  axis  Actually  only  the  positive 

or  the  negative  spectrum  need  be  used,  for  the  omission  of  either  merely 
causes  a light  loss  by  a factor  of  two  in  the  intensity  of  the  correlation  function 
in  the  output  plane  This  result  is  demonstrated  in  Appendix  If 

In  general,  the  signal  spectrum  is  complex,  and,  therefore,  both 
amplitude  (opacity)  and  phase  (thickness)  recordings  must  be  included  in  the 
optical  version  of  the  reference  function  inserted  in  plane  (Fig  3) 

However,  there  is  a large  class  of  signals  for  which  no  phase  variation  is 
required  These  include  all  signals  whose  Fourier  transforms  are  purely  real 
or  purely  imaginary,  corresponding  to  signals  which  are  purely  even  or  purely 
odd  functions  of  time  The  use  of  such  signals  considerably  simplifies  the 
construction  of  an  optical  reference,  for  only  amplitude  recording  is  then 
necessary 

When  the  inverse  transform  is  taken  by  lens  L2  of  Fig  3,  a sum  of 
superimposed  correlation  functions,  one  for  each  light  frequency,  will  appear 
in  plane  P3  A more  exact  description  would  note  that  the  optic  »xis  for 
lens  L2  in  Fig  3 cannot  be  colinear  with  the  diffracted  optic  axes  for  ail 
wavelengths  This  lack  of  colinearity  has  two  effects  upon  the  correlation 
function  in  plane  P3: 

1)  A slight  displacement  of  diffracted  axes,  relative  to  the 

axis  for  lens  L2  and  other  components  to  the  right  of  L2  in  Fig  3,  exists 
at  all  light  wavelengths  except  for  a particular  matched  v velengih  X This 
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small  displacement  in  a frequam  y shift  which  cun  hr  ignored  however  mni  r 

a (llapUcement  »X,*>^(X'  in  the  positive  ^or  negative!  frequency  plane  produces 

* multiplying  phase  factor  exp  [- j ( X;  x ^ | in  the  next  spatial  plane  p.i 

where  the  correlation  function  is  imaged  Since  light  intensity  rather  than 

amplitude  is  sensed  by  an  output  film  or  photodetector  the  absolute  nquarr 

of  the  correlation  coefficient  is  detected  and  phase  factors  of  the  t\pe 

exp  ( , K!  x ^ ] are  not  apparent 

Z)  in  general,  the  plane  P«J^  for  w given  light  frequency  is  slightly 

inclined  (not  exactly  perpendicular)  to  the  optic  axis  of  lens  L 2 Thus  the 

plane  P2  o which  is  actually  perpendicular  to  the  L2  axis  and  in  which  the 

reference  function  is  inserted  has  a light  amplitude  distribution  which  for  atb» 

trary  X ta  not  quite  the  Fourier  transform  of  that  existing  in  lane  FI  but 

rather  Is  one  which  bv  reason  of  the  *tlted  wavefront  P2^  can  he  described  by 

assigning  to  the  Fourier  transform  « phase  shift  linear  with  displacement  in  t'.v 

F2vo  plane  Thua  the  light  amplitude  distribution  transmitted  through  P2  o ia 

Jka)u.  J*:X!o  K 

«il  the  form  Sfw^'  R*(w^)  e x where  e i«  the  phase  factor 

and  the  const  .nt  k(X)  is  a measure  i><  the  angle  between  P2^  and 

The  action  of  lens  L2  in  performing  the  inverse  Fourier  transform  converts 

this  phase  factor  to  a time  displacement  ktX)  in  plane  * and  this  time 

displacement  is  dependent  upon  \ Thus  the  time  skis  in  Pf  does  not 

have  the  same  origin  for  all  X and  the  output  plane  at  P4  must  be  prop 

erly  calibrated  with  regard  to  X 


P M Woodward  Probabiltt y and  rdormatinn  Theory  with  Applications 
to  Radar  (McOnw  Hill  Rook  C©  Inc  . New  York  iqx ) t «ble  \ 20 


The  light  intensity  available  at  plane  P2^  ♦'hen  ilium  nated 
fr<  m the  on-axis  point  of  plane  PI  is  given  by  the  usual  expression  for 
the  l)ght  intensity  distribution  from  a plane  wavefront  mpngng  on  a grating 


^ 2 jiin^p  sin^Ny 


sin^Y 
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where 


C 


P = 
0 = 

a - 

Y - 

d - 


it  a sin  0 
\ 

diffraction  angle 

grating  slit  width 
it  d sin  6 


grating  slit  separation,  and 
N - number  of  slits  in  the  grating 

By  making  the  "lit  width  a large  (e  g equal  to  one  h *lf  the  slit  eparation 

sin^e 

d it  is  possible  to  narrow  the  single  slit  modulating  envelope  

P 

so  that  most  of  the  light  is  concentrated  in  the  n - 0 ± 1 orders  There 
is  also  a requirement  that  the  various  diffraction  orders  do  not  overlap  at 
plane  P2  and  thereby  confuse  the  re«ulting  image  This  requirement  is 
met  by  using  long  focal  length  lense*  and  small  slit  separation  d on  the 
grating  to  provide  s large  diffraction  angle  0 The  joint  requirements 
of  light  intensity  distribution  and  image  separation  can  be  achieved  by  making 
d 2a  with  d small  It  is  also  possible  on  ruled  gratings  t-  t»per  the 
bottom  of  each  cut  so  as  to  concentrate  up  to  90  percent  o<  the  transmitted 
light  intensity  into  the  first  {e  g n - 1 / order  at  the  expense  o the  central 
(n  ~ 0)  and  higher  order  images  Phase  gratings  in  which  all  parts  o(  the 
grating  are  equally  transparent  but  difrer  in  optical  thickness  net  onlv  preserve 
that  light  normally  lost  in  the  opaque  regions  of  convent  oral  gratings  but  can 
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aide  bo  constructed  to  concentrate  light  in  the  first  order  images  at  the 
expense  of  the  central  image 

The  diffraction  pattern  from  grating  G1  has  secondary  maxima  at 
angles  close  to  the  primary  maxima  and  these  w ill  add  incoherent  illumination 
to  the  primary  maximum  diffraction  image  in  plane  P2^  It  can  be  shown 
however  that  the  intensity  of  a particula  secondary  maximum  relative  to 
that  oi  the  primary  maximum  is  approximately  4 (2n  -t  l)^  where 

n t6  the  number  of  the  secondary  maximum  c ounting  from  the  nearest 
prima-y  maximum  The  total  intensity  received  from  Mil  secondary  maxima 
within  an  angular  width  ±0  c>/2  , the  assumed  angular  width  of  the  image, 

to  either  side  of  a primary  maximum  is  approximately 


8 

~T 


(N-2)/2 

n?i 


1 


20 


(2n-t  It 


(12. 


However  the  angular  width  of  the  principal  maximum  is  twice  that  of  the 
secondaries  Therefore  the  relative  noncoherent  illuminati  n supplied  by 
the  secondary  maxima  is  approximately  one  half  of  Eq  (12)  or  10  percent 
of  the  intensity  in  the  principal  maximum  In  most  case**  the  10  percent 
increase  in  n >i»e  level  will  not  significantly  affect  the  detectability  of  * 
signal  which  usually  is  recorded  with  additive  noise  o'  much  larg-r  pro- 
portions 

it  might  seem  that  the  angular  width  of  the  n t 1 or  - 1 principal 
maximum  would  also  limit  system  performance  by  degrading  the  inherent 
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angular  revolution  However,  it  is  easily  shown’  that  this  is  an  aperture 
effect,  independent  of  the  grating  constants  and  so  lonR  as  the  aperture  in 
which  the  grating  is  placed  is  the  same  sire  as  other  apertures  of  the 
system  the  angular  resolution  will  be  the  same  throughout 

The  prism  preceding  the  output  plane  P4  of  Fig  3 separates  the 
various  li^ht  colors,  or  Doppler  channel*  into  a continuous  spectrum  spread 
vertically  This  vertical  spread  need  not  produce  any  overlap  of  closely 


The  angular  location  of  the  first  minimum  to  either  side  of  the  n ♦ 1 
principal  maximum  i*  given  by  ’ 


sin  0 


min 


jN  ± 1'  X_ 
Nd 


while  that  of  the  principal  maximum  is 


sin  0 


max 


For  large  N the  half  width  ; 0 - 6 . of  the  principal  maximum 

is  obtained  from  rnax  01,0 


■in  ® sin  0 

max  min 


d (gin  0 ) 

max 


- X / Nd 


cos  0 d0 

max  max 


d9  * X/  Nd  cos  0 B X A 

max  max 


where  A = total  aperture  This  is  the  approximate  expression  for  the 
angular  resolution  of  a circular  aperture  of  diameter  A 
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spaced  signed  channels  from  the  signal  film  ‘plane  FI),  for  optical  magni 
fication  between  plane3  P3  and  P4  can  produce  many  more  resolvable 
channels  on  P4  than  existed  in  PI  by  enlarging  the  format  Moreover 
the  vertical  or  y^  , spacing  of  channels  in  PI  is  dictated  by  the  overall 
resolution  of  the  correlator  as  * whole  (film  plus  lens  components  acting  as 
cascaded  filters),  whereas  the  channel  spacing  in  plane  P4  need  only  exceed 
the  resolution  distance  on  the  film  alone  It  should  be  noted  that  the  optical 
transformation  between  planes  P3  and  P4  is  a true  imaging  process  and 
not  a Fourier  transform  relationship  Put  in  different  terms  the  power  of 
lens  L3  is  such  that  both  the  Fourier  transform  and  its  inverse  are  accom- 
plished between  P3  and  P4  , 

The  prism  can  also  be  oriented,  as  shown  in  Fig  6,  so  as  to  produce 
a horizontal,  or  x , spread  of  the  spectrum  To  avoid  confusion  between  the 
time  axis  and  Doppler  channels,  it  would  then  be  necessary  to  insert  a slit  in 
plane  P3  and  use  the  x axis  in  plane  p4  for  the  Doppler  X)  as  veil  as 
the  time  <x0'  coordinate  as  shown  in  Fig  6 In  this  configuration  the 
correlation  function  <J>  (x^>  is  not  simultaneously  observable  o\er  a wide  x 
range  as  it  is  in  the  configuration  of  Fig  2 instead,  only  <t>  x 1 for  that 
Xo  v*lue  corresponding  to  the  slit  position  is  sampled  at  one  inst  »nt  of  time 
With  the  signal  film  continuously  traversing  plane  PI  in  real  time,  the  out- 
put through  the  slit  at  plane  P3  (or  its  Doppler  separated  version  at  P4; 
provides  < time  evolving  history  of  4>  lx„  To  avoid  information  smear  on 
the  output  film,  the  Doppler  spread  l\  \ ) must  be  sample'"  at  P4  in 
time  increments  A xq  . and  the  Doppler  determinations  must  be  multiplexed 
with  the  real  time  x^  variation  This  leads  naturalky  to  some  sort  of 


frame  camera  for  recording  information  appearing  at  P4  but  the  resulting 
multiplexing  of  Doppler  and  time,  or  , information  may  not  be  particu- 
larly convenient  for  subsequent  data  reduction  Moreover,  a much  shorter 
exposure  time,  equal  to  £ is  available  in  the  configuration  of  Fig  6 

than  in  that  of  Fig  3,  where  a plot  of  ♦(xq)  for  a range  of  xq  values  is 
available  at  P3  (or  P4)  and  exposes  a film  moving  slowly  through  the  P3 
(or  P4)  plane  synchronously  with  the  signal  film  motion  through  PI 

For  some  purposes  the  prism  can  be  eliminated  entirely  and  the  various 
Doppler  channels  left  unresolved  in  plane  P4  The  correlator  then  becomes 
Doppler  invariant  regardless  o«‘  whether  ir  not  the  original  s gn»l  function 
was  chosen  to  be  Doppler  invariant  Th  e >«n  be  a very  c men  ent  mode  of 
operation  when  one  wishes  not  only  to  make  no  use  of  Doppler  nft  mution  but 
*lso  needs  to  insure  that  Doppler  distorticn3  of  a target  reflected  signal  do 
not  cause  the  original  signal  to  decorrelate  Ot  course,  a particular  constant 
velocity  target  motion  will  then  give  rise  to  a received  signal  which  correlates 
at  only  one  light  frequency  and  illumination  in  the  output  plane  from  all  other 
parts  of  the  light  spectrum  will  then  appear  as  noise  superimposed  on  the 
desired  correlation  function  Therefore  if  the  Doppler  channels  are  left 
unresolved  in  P4  , the  signal  to  noise  ratio  in  plane  PI  must  be  adequate 
to  insure  a detection  at  P4  in  the  presence  of  noise  from  other  light  frequencies 
Since  the  distortion  of  the  signal  timebase  by  a scale  factor  due  to 
moving  target  reflections,  or  the  distortion  of  its  frequency  components  by 
the  inverse  of  this  scale  factor,  is  normally  a small  percentage  distortion, 
the  range  of  light  frequencies  required  for  Doppler  compensation  ie  a small 
percentage  of  the  center  frequency  .see  Eq  [ 1 ] > Therefore,  no  compensation 
need  be  made  for  the  spectral  response  of  the  light  source  or  the  film  in  the 
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output  plane  P4  over  the  narrow  apectral  range  actual' y ueed  However, 
the  tame  usually  cannot  be  said  for  the  band-pass  interference  filter  FI 
(Fig  3),  whose  frequency  response  cuvve  has  skirts  that  do  not  fall  to  r.ero 
infinitely  fast  at  the  edges  of  the  pass  region  To  the  extent  that  compen 
sation  must  be  made  the  light  intensities  in  plane  P4  m»y  bo  interpreted 
after  weighting  them  with  the  inverse  frequency  response  of  the  filter  FI 

The  relative  merits  of  the  panchromatic  principle  can  be  summarised 
by  reviewing  the  advantages  and  disadvantages  implied  in  the  foregoing  The 
advantages  relative  to  other  optical  correlation  techniques  appear  to  include 
the  following: 

1)  A single  reference  function  for  correlation  against 
all  Doppler  distortions  of  the  signal 

2)  Simultaneous  correlation  of  all  Doppler  distorted 
signals,  rather  than  a sequential  Doppler  search 
multiplexed  with  real  time  (or  range) 

3)  No  mechanical  motions  required  either  for  Doppler 
search  or  heterodyning 

4)  Longer  time  available  to  expose  film  in  the  output 
plane  P4  , by  reason  of  2)  above. 

**)  Simplicity  of  parts  relative  to  the  number  and  com 
plexity  of  operations  accommodated  (heterodyning 
Doppler  compensation,,  correlation) 

The  principal  disadvantage  of  the  panchromatic  configuration  results 
from  its  inefficient  use  of  light  flux  available  at  the  source  Roughly  90  per- 
cent of  the  initial  light  flux  is  loat  by  retaining  only  one  order  n value)  In  the 
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diffraction  image  produced  by  simple  gratings  Where  light  source  intensity 
is  not  a limitation,  or  where  advantage  4)  , above,  compensates,  or  when 
special  gratings  are  used  to  concentrate  light  in  the  n = - 1 order,  this  is 
obviously  not  a problem  In  any  event,  the  potential  capability  of  the  pan- 
chromatic configuration  to  perform  simultaneous  Doppler  compensation  and 
correlation  operations  on  heterodyned  signals  appears  intriguing  and  theo 
retically  promising 
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APPEND  X X 


Ihf  cc-mplex  light  amplitude  v^ct  r &2  <X2.V?)  at  point  x^  . y 
of  plane  P2  (ser  Fig  7 } result*  from  summing  the  contributions  Ej'x]  . Yj) 
from  all  points  'Xj  , y^  of  plana  PI  Neglecting  the  obliquity  factor 
and  the  time  dependent  factor  thin  vector  sum  can  be  written  ^ 


^2  'KZ‘fZ)  ~ fjE  1(X1  yi;  e 


j-rr(xi  vi  x2  V 


dxv  dy,  , (13) 


where  rlXj  ,yj  ’ yZ  *8  -he  cptical  path  length  between  points  x^  Yj) 
in  PI  and  x^  y in  P2  When  a dif  rsction  grating  is  present,  Eq  <13 
can  be  used  to  calculate  the  amplitude  vector  in  any  order  of  the  diffraction 
image  providing  r represents  the  distance  to  the  image  point  in  that  order 
and  the  various  orders  are  non  o^  erlapping  Figure  chews  the  geometry 
involved  for  a particular  wavelength  and  the  n J first  order  diffracted 
image  Xt  is  apparent  from  the  figure  th*t  the  optical  distance  r x,  , from 
any  point  of  the  tilted  plane  PI'  to  is 


o *2 


tj  •i  f 


, = g COB  i r — , 

2 ° cob v9  0 t 


where  is  the  incident  ray  and  r2  the  diffracted  ray  to  the  first  order 
diffracted  image  The  condition  for  the  first  order  principal  maximum  of 


the  diffraction  im?ge  ia 


d (sin  8 sin  i)  X 


sin  i r.  gin  8 - gin  8 


fr-’T-cry 

H! ' m 
t :du(84) 
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where  a n 0 — j and  0^  is  the  angular  position  of  the  princjp*JL 

maximum  v»h»n  ».  o Substituting  for  cos  Eq  (14)  becomes 


ro<*2) 


(3in0  einG^  ^ tan^(0  0 J 


Expanding  the  trigonometric  functions  and  retaining  powers  of  the  angles 
below  the  fourth 
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Any  other  point  o;  plane  PI  will  be  closer  to  or  further  from  than 


plane  PI'  by  a distance  - x ^ sin  i ) 
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whare  the  second  term  on  tie  right  aide  involving  third  order  differences  in 
0 rnd  0 Rnd  their  products,  is  neglected 

To  justify  this  neglect  it  must  fc>e  shown  that  the  value  of  this  term  is 
of  the  same  order  or  smaller  than  the  angular  resolution  due  to  the  aperture 
ie»  the  angular  width  of  a principal  maximum  in  the  diffraction  image t which 
is  given  by  X.  / Nd  t is  not  hard  to  choose  a set  of  practical  values  for  N,  d 
and  f which  meet  this  condition  while  still  avoiding  an  overlap  of  adjacent  orders 
in  the  diffraction  image 

Adding  Eq  v 1 6)  to  Eq.  (15' 


x 1 x2  Xh 


r ixj  4 g + f + *f  - g i io  -0ol2  • -*  - f - 


Allowing  g - f and  neglecting  the  constant  teems  g and  f the 

remaining  term  is  just  - k.  ix  ~x  } f For  two  dimensions  the  optical 

i w n 

distance  r^x^y^  y^j  , apart  from  constant  terms,  becomes 
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Equation  13  can  new  be  written 
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- 

*V  = ’ xF  v2 


Since  x?  xhj  and  y,  , and  therefore  wj;  anc  . are  measured  from 

the  di  > *'*ed  optic  ?x;,8  whose  location  & a function  of  X t^ey  will  have 

di^e  rent  values  when  describing  a particular  paint  in  plane  P2  o forth? 

\ 

various  lig.b  wavelenptha  X~  \ " 
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APPEND  X 


The  Fourier  t ■Hnsf'orm  o:  the  delayed  r-s  :eived  sign  1 fft»  ' can  be 


written 
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additive  no  se  positive  frequency 
spectrum  i i the  received  function  fit}. 

positive  fr<  quen-y  spectrum  of  the 
signal  s (It 

constant  <<  1 representing  signal 
amplitude  r elative  to  the  noise, 

positive  frequency  spectrum  of  the 
teceiv  ed  function  f(t),  ’end 

positive  frt  quency  spectrum  nf  the 
delayed  ret  eived  function  f c 1 4> 


The  positive  frequency  spectrum  obiaine  1 from  Eq  1 8 r after  hetero- 
dyning by  fin  amount  2 is 
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X 

and  this  quantity  is  multiplied  by  S ;w  SI  in  ? mat  hed  filter  receiver, such 
an  the  panchromatic  correlator  to  form 
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However  it  both  poaitive  and  nemfiie  frequencies  h d bpi>n  retained 
m the  matched  liter  operation  the  correlation  function  mould  h ■ v ,•  been 


♦ ‘ 2„  { Sv 


But  t»  • an>  e »1  signal  e(t) 


2 i - u»  , U 

e ' 


I »*>  t K 

e • dui  *-  S 
2 n J 

UO 


11 


’ J ll!  jj 

f 0*0 


S+(U)40  ) S ( - u)  11 


s^wffl);2  S ( W n ) 1 


Thu  a 


y f |S  i u ! 11  | ^ J i ^ f ii  V'  ♦ j<*>t  e J(«  11 
0 + 


K /'  2 

I S juii  Q ) [ cn  u;t»',  coa  il  ain  .it  >ln  u ] d o 


* 2 f X lt+  i rof  11  ' Y ( t ' i ' ain  W t 


♦ 2^X(t+v)j2  ♦ 21  Y ft  p ' i 2 { f X t )f  | V (t  '*}< 


« 2 n 


x tt-t  | Y {t*  t ) ain  211 


Because  the  heterodyning  frequency  shift  ft  is  generally  larger  than  the 
signal  bandwidth  accepted  by  the  optical  correlator,  such  terms  as  cos  2ft  r 
and  sin  2ft  t average  to  rero  within  the  space  covered  by  one  resolvable  spot 
on  the  output  film  (or  photodetector)  Then  the  light  intensity  which  is  actually 


detected  becomes 


just  twice  the  value  of  «J>  (t+r)  j Thus  the  net  effect  of  operating  only 

on  the  positive  frequency  spectrum  of  the  received  function  f(t)  is  one  of 
reducing  the  available  light  level  at  the  output  by  a factor  of  two  For  mechanical 
convenience  the  panchromatic  correlator  is  designed  to  accept  only  the  positive 

i i2 

frequency  spectrum,  thereby  generating  A (t+T)  | at  its  output 
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